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O ■ Abstract 

\ Black holes came into existence together with the universe through the 

quantum process of pair creation in the inflationary era. We present the 
q-i instantons responsible for this process and calculate the pair creation rate 

5^ \ from the no boundary proposal for the wave function of the universe. We find 

that this proposal leads to physically sensible results, which fit in with other 
descriptions of pair creation, while the tunnelling proposal makes unphysical 
^ \ predictions. 

' We then describe how the pair created black holes evolve during inflation. 

In the classical solution, they grow with the horizon scale during the slow 
roll-down of the inflaton field; this is shown to correspond to the flux of field 
energy across the horizon according to the First Law of black hole mechanics. 
When quantum effects are taken into account, however, it is found that most 
black holes evaporate before the end of inflation. Finally, we consider the 
pair creation of magnetically charged black holes, which cannot evaporate. In 
standard Einstein-Maxwell theory we find that their number in the presently 
observable universe is exponentially small. We speculate how this conclusion 
may change if dilatonic theories are applied. 
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1 Introduction 



It is generally assumed that the universe began with a period of exponential ex- 
pansion called inflation. This era is characterised by the presence of an effective 
cosmological constant A e ff due to the vacuum energy of a scalar field <fi. In generic 
models of chaotic inflation Q , the effective cosmological constant typically starts 
out large and then decreases slowly until inflation ends when A e g- m 0. Correspond- 
ingly these models predict cosmic density perturbations to be proportional to the 
logarithm of the scale. On scales up to the current Hubble radius this agrees 

well with observations of near scale invariance. However, on much larger length 
scales of order H~^ w exp(10 5 ), perturbations are predicted to be on the order of one. 
Of course, this means that the perturbational treatment breaks down; but it is an 
indication that black holes may be created, and thus warrants further investigation. 

An attempt to interpret this behaviour was made by Linde || |J . He noted that 
in the early stages of inflation, when the strong density perturbations originate, the 
quantum fluctuations of the inflaton field are much larger than its classical decrease 
per Hubble time. He concluded that therefore there would always be regions of the 
inflationary universe where the field would grow, and so inflation would never end 
globally ("eternal inflation"). However, this approach only allows for fluctuations 
of the field. One should also consider fluctuations which change the topology of 
space-time. This topology change corresponds to the formation of a pair of black 
holes. The pair creation rate can be calculated using instanton methods, which are 
well suited to this non-perturbative problem. 

One usually thinks of black holes forming through gravitational collapse, and 
so the inflationary era may seem an unlikely place to look for black holes, since 
matter will be hurled apart by the rapid cosmological expansion. However, there 
are good reasons to expect black holes to form through the quantum process of pair 
creation. We have already pointed out the presence of large quantum fluctuations 
during inflation. They lead to strong density perturbations and thus potentially 
to spontaneous black hole formation. But secondly, and more fundamentally, it is 
clear that in order to pair create any object, there must be a force present which 
pulls the pair apart. In the case of a virtual electron-positron pair, for example, 
the particles can only become real if they are pulled apart by an external electric 
field. Otherwise they would just fall back together and annihilate. The same holds 
for black holes: examples in the literature include their pair creation on a cosmic 
string []|], where they are pulled apart by the string tension; or the pair creation of 
magnetically charged black holes on the background of Melvin's universe || , where 
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the magnetic field prevents them from recollapsing. In our case, the black holes will 
be separated by the rapid cosmological expansion due to the effective cosmological 
constant. So we see that this expansion, which we naively expected to prevent black 
holes from forming, actually provides just the background needed for their quantum 
pair creation. 

Since inflation has ended, during the radiation and matter dominated eras until 
the present time, the effective cosmological constant was nearly zero. Thus the 
only time when black hole pair creation was possible in our universe was during the 
inflationary era, when A e g- was large. Moreover, these black holes are unique since 
they can be so small that quantum effects on their evolution are important. Such 
tiny black holes could not form from the gravitational collapse of normal baryonic 
matter, because degeneracy pressure will support white dwarfs or neutron stars 
below the Chandrasekhar limiting mass. 

In the standard semi-classical treatment of pair creation, one finds two instan- 
tons: one for the background, and one for the objects to be created on the back- 
ground. From the instanton actions ibg and 7 bj one calculates the pair creation rate 

r: 

r = exp [- (/ obj - / bg )] , (1.1) 

where we neglect a prefactor. This prescription has been very successfully used by 
a number of authors recently |S], [II], |TT], [12], for the pair creation of black 

holes on various backgrounds. 

In this paper, however, we will obtain the pair creation rate through a somewhat 
more fundamental, but equivalent procedure: since we have a cosmological back- 



ground, we can use the Hartle-Hawking no boundary proposal |14| for the wave 
function of the universe. We will describe the creation of an inflationary universe by 
a de Sitter type gravitational instanton, which has the topology of a four-sphere, S* 4 . 
In this picture, the universe starts out with the spatial size of one Hubble volume. 
After one Hubble time, its spatial volume will have increased by a factor of e 3 ~ 20. 
However, by the de Sitter no hair theorem, we can regard each of these 20 Hubble 
volumes as having been nucleated independently through gravitational instantons. 
With this interpretation, we are allowing for black hole pair creation, since some 
of the new Hubble volumes might have been created through a different type of 
instanton that has the topology S 2 x S 2 and thus represents a pair of black holes in 
de Sitter space ||15|| . Using the framework of the no boundary proposal (reviewed in 
Sec. H), one can assign probability measures to both instanton types. One can then 
estimate the fraction of inflationary Hubble volumes containing a pair of black holes 



3 



by the fraction T of the two probability measures. This is equivalent to saying that 
T is the pair creation rate of black holes on a de Sitter background. We will thus 
reproduce Eq. (|1 . 1|) . 

In Sec. |3.1| we follow this procedure using a simplified model of inflation, with a 
fixed cosmological constant, before going to a more realistic model in Sec. |3.2| . In 
Sec. |3.3j we show that the usual description of pair creation arises naturally from the 



no boundary proposal, and Eq. (|1 . 1| ) is recovered. We find that Planck size black 



holes can be created in abundance in the early stages of inflation. Larger black 
holes, which would form near the end of inflation, are exponentially suppressed. 
The tunnelling proposal fIB |, on the other hand, predicts a catastrophic instability 
of de Sitter space and is unable to reproduce Eq. ( p..l|) . 

We then investigate the evolution of black holes in an inflationary universe. 
In Sec. f| their classical growth is shown to correspond to energy-momentum flux 
across the black hole horizon. Taking quantum effects into account, we find in 
Sec. |5| that the number of neutral black holes that survive into the radiation era is 
exponentially small. On the other hand, black holes with a magnetic charge can 
also be pair created during inflation. They cannot decay, because magnetic charge 
is topologically conserved. Thus they survive and should still be around today. In 
Sec. |6], however, we show that such black holes would be too rare to be found in 
the observable universe. We summarise our results in Sec. |7|. We use units in which 
mp = h = c = k = l. 



2 No Boundary Proposal 

We shall give a brief review; more comprehensive treatments can be found else- 



where [pL7| . According to the no boundary proposal, the quantum state of the 
universe is defined by path integrals over Euclidean metrics g^ on compact mani- 
folds M. From this it follows that the probability of finding a three-metric on a 
spacelike surface £ is given by a path integral over all g^ u on M that agree with 
on E. If the spacetime is simply connected (which we shall assume), the surface £ 
will divide M into two parts, M + and M_. One can then factorise the probability of 
finding into a product of two wave functions, and ^> -. _) is given by a 

path integral over all metrics g^ v on the half-manifold M + (M_) which agree with 
hij on the boundary E. In most situations equals We shall therefore drop 
the suffixes and refer to \1/ as the wave function of the universe. Under inclusion of 
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matter fields, one arrives at the following prescription: 

*e] = / $) exp [-I{ 9lxv , $)] , (2.1) 

where (/iy , $s) are the 3- metric and matter field on a spacelike boundary £ and the 
path integral is taken over all compact Euclidean four geometries g^ that have S 
as their only boundary and matter field configurations $ that are regular on them; 
I(<7 M „, $) is their action. The gravitational part of the action is given by 

I E = — / d 4 x g 1/2 (R - 2A) - — / Sx h l/2 K } (2.2) 

167T JM+ 87T JS 

where R is the Ricci-scalar, A is the cosmo logical constant, and K is the trace of 
Kij, the second fundamental form of the boundary £ in the metric g. 

The wave function \l/ depends on the three-metric hij and on the matter fields $ 
on E. It does not, however, depend on time explicitly, because there is no invariant 
meaning to time in cosmology. Its independence of time is expressed by the fact 
that it obeys the Wheeler-DeWitt equation. We shall not try to solve the Wheeler- 
DeWitt equation directly, but we shall estimate $ from a saddle point approximation 
to the path integral. 

We give here only a brief summary of this semi-classical method; the procedure 
will become clear when we follow it through in the following section. We are inter- 
ested in two types of inflationary universes: one with a pair of black holes, and one 
without. They are characterised by spacelike sections of different topology. For each 
of these two universes, we have to find a classical Euclidean solution to the Einstein 
equations (an instanton), which can be analytically continued to match a boundary 
£ of the appropriate topology. We then calculate the Euclidean actions I of the two 
types of solutions. Semiclassically, it follows from Eq. ( |2.1|) that the wave function 
is given by 

* = exp(-J), (2.3) 

where we neglect a prefactor. We can thus assign a probability measure to each type 
of universe: 

P= |^| 2 = exp(-2/ Re ) , (2.4) 

where the superscript 'Re' denotes the real part. As explained in the introduction, 
the ratio of the two probability measures gives the rate of black hole pair creation 
on an inflationary background, T. 
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3 Creation of Neutral Black Holes 



The solutions presented in this section are discussed much more rigorously in an 
earlier paper [|lSj . We shall assume spherical symmetry. Before we introduce a more 



realistic inflationary model, it is helpful to consider a simpler situation with a fixed 
positive cosmological constant A but no matter fields. We can then generalise quite 
easily to the case where an effective cosmological "constant" arises from a scalar 
field. 

3.1 Fixed Cosmological Constant 
3.1.1 The de Sitter solution 

First we consider the case without black holes: a homogeneous isotropic universe. 
Since A > 0, its spacelike sections will simply be round three-spheres. The wave 
function is given by a path integral over all metrics on a four-manifold M + bounded 
by a round three-sphere £ of radius as- The corresponding saddle point solution 
is the de Sitter space-time. Its Euclidean metric is that of a round four-sphere of 
radius J 3/A: 

ds 2 = dr 2 + a(r) 2 dnl (3.1) 

where r is Euclidean time, dVt\ is the metric on the round three-sphere of unit 
radius, and 

a(T) = /f Sin /f r - (3 - 2) 
For a s = 0, the saddle point metric will only be a single point. For < a s < 



3/A it will be part of the Euclidean four-sphere , and when as = J 3/A, the saddle 



point metric will be half the four-sphere. When as > y 3/A there will be no real 
Euclidean metric which is a solution of the field equations with the given boundary 
conditions. However, we can regard Eq. ( |3.2| ) as a function on the complex r-plane. 
On a line parallel to the imaginary r-axis defined by r Re = J f; § , we have 



•1 JL 

A 2 



'l coslJV*. (3.3) 



This describes a Lorentzian de Sitter hyperboloid, with r Im serving as a Lorentzian 
time variable. One can thus construct a complex solution, which is the analytical 
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continuation of the Euclidean four-sphere metric. It is obtained by choosing a 
contour in the complex r-plane from to r Rc = \J\\ and then parallel to the 
imaginary r-axis. One can regard this complex solution as half the Euclidean four- 
sphere joined to half of the Lorentzian de Sitter hyperboloid (Fig. 1). 




Figure 1: The creation of a de Sitter universe. The lower region is half of a Euclidean 
four-sphere, embedded in five-dimensional Euclidean flat space. The upper region 
is a Lorentzian four-hyperboloid, embedded in five-dimensional Minkowski space. 

The Lorentzian part of the metric will contribute a purely imaginary term to the 
action. This will affect the phase of the wave function but not its amplitude. The 
real part of the action of this complex saddle point metric will be the action of the 
half Euclidean four-sphere: 

ifte Sitter = ~2A' fi-ty 
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Thus the magnitude of the wave function will still be e 3?r//2A , corresponding to the 
probability measure 

Pdc sitter = exp (-^J . (3.5) 



A 

3.1.2 The Schwarzschild-de Sitter solution 

We turn to the case of a universe containing a pair of black holes. Now the cross 
sections £ have topology S 2 x S 1 . Generally, the radius of the S 2 varies along the S 1 . 
This corresponds to the fact that the radius of a black hole immersed in de Sitter 
space can have any value between zero and the radius of the cosmological horizon. 
The minimal two-sphere corresponds to the black hole horizon, the maximal two- 
sphere to the cosmological horizon. The saddle point solution corresponding to this 
topology is the Schwarzschild-de Sitter universe. However, the Euclidean section of 
this spacetime typically has a conical singularity at one of its two horizons and thus 
does not represent a regular instanton. This is discussed in detail in the Appendix. 
There we show that the only regular Euclidean solution is the degenerate case where 
the black hole has the maximum possible size. It is also known as the Nariai solution 
and given by the topological product of two round two-spheres: 

ds 2 = dr 2 + a(rfdx 2 + b(r) 2 dn 2 2 , (3.6) 

where x is identified with period 2ir, dil 2 . = d6 2 + sin 2 6> dip 2 , and 

a(r) = ]j^smy/Xr, b(r) = = const. (3.7) 

In this case the radius b of the S 2 is constant in the S 1 direction. The black hole 
and the cosmological horizon have equal radius 1 / y/X and no conical singularities are 
present. Thus, by requiring the smoothness of the Euclidean solution, the instanton 
approach not only tells us about probability measures, but also about the size of the 
black hole. There will be no saddle point solution unless we specify 6 S = I /y/X. We 
are then only free to choose the radius as of the one-sphere on E. For this variable, 
the situation is similar to the de Sitter case. There will be real Euclidean saddle 
point metrics on M + for as < 1/ y/X. For larger there will again be no Euclidean 
saddle point, but we find that 

aWI^/r, = \ljCOshy/Xr lm . (3.8) 



A 2 
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This corresponds to the Lorentzian section of the degenerate Schwarzschild-de Sitter 
spacetime, in which the S 1 expands rapidly, while the two-sphere (and therefore the 
black hole radius) remains constant. Again we can construct a complex saddle point, 
which can be regarded as half a Euclidean S 2 x S 2 joined to half of the Lorentzian 
solution. The real part of the action will be the action of the half of a Euclidean 
S 2 x S 2 : 

's R dS = ~- (3-9) 
The corresponding probability measure is 

Psds = exp . (3.10) 



We divide this by the probability measure ( |3.5|) for a universe without black holes 
to obtain the pair creation rate of black holes in de Sitter space: 

-Psds / TT \ 

ex P _ T ■ ( 3 - n ) 



dc Sitter 



A 



Thus the probability for pair creation is very low, unless A is close to the Planck 
value, A = 1. 



3.2 Effective Cosmological Constant 

Of course the real universe does not have a cosmological constant of order the Planck 
value. However, in inflationary cosmology it is assumed that the universe starts out 
with a very large effective cosmological constant, which arises from the potential V 
of a scalar field 0. The exact form of the potential is not critical. So for simplicity we 
chose V to be the potential of a field with mass m, but the results would be similar 
for a A0 4 potential. To account for the observed fluctuations in the microwave 
background ||19|| , m has to be on the order of 10~ 5 to 10 -6 [pOj] . The wave function 
\1/ will now depend on the three-metric and the value of on S. By a gauge 
choice one can take to be constant on E, and we shall do so for simplicity. For 
> 1 the potential acts like an effective cosmological constant 

A cfr (0) = 8ttV(0). (3.12) 

One proceeds in complete analogy to the fixed cosmological constant case. For 
small three-geometries and > 1, there will be real Euclidean metrics on M+, with 
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almost constant. If the three-geometries are rather larger, there will again not 
be any real Euclidean saddle point metrics. There will however be complex saddle 
points. These can again be regarded as a Euclidean solution joined to a Lorentzian 
solution, although neither the Euclidean nor Lorentzian metrics will be exactly real. 
Apart from this subtlety, which is dealt with in Ref. |18[], the saddle point solutions 
are similar to those for a fixed cosmological constant, with the time-dependent A e g 
replacing A. The radius of the pair created black holes will now be given by 1/ y/K^g. 
As before, the magnitude of the wave function comes from the real part of the action, 
which is determined by the Euclidean part of the metric. This real part will be 



rRe 
J S 3 



in the case without black holes, and 



3tt 



2A cff (0 o ) 



(3.13) 



t Re 

1 S 2 xS 1 



71 



A, 



eff^Oj 



in the case with a black hole pair. Here 0o is the value of 
region. Thus the pair creation rate is given by 



(3.14) 



in the initial Euclidean 



P S 2 



xS 1 



P s s 



cxp 



7T 



A e ff(0o) 



(3.15) 



3.3 Discussion 

Let us interpret this result. Since < A e g < 1, we get T < 1 and so black hole pair 
creation is suppressed. In the early stages of inflation, when A e g « 1, the suppression 
is week, and black holes will be plentifully produced. However, those black holes 
will be very small, with a mass on the order of the Planck mass. Larger black holes, 
corresponding to lower values of A e g- at later stages of inflation, are exponentially 
suppressed. We shall see in the following two sections that the small black holes 
typically evaporate immediately, while sufficiently large ones grow with the horizon 
and survive long after inflation ends (that is, long in early universe terms). 

We now understand how the standard prescription for pair creation, Eq. Q1.1J ), 
arises from this proposal: by Eq. (|2.4j ), 



S 2 xS 1 

P 5 3 



exp 



(21 



Re 

S 2 xS 1 



21}! 



(3.16) 
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where J Re denotes the real part of the Euclidean action of a complex saddle point 
solution. But we have seen that this is equal to half of the action of the complete 
Euclidean solution. Thus J bj = 2lfi xSl and = 21 g§, and we recover Eq. (|1.1|). 

The prescription for the wave function of the universe has long been one of the 
central, and arguably one of the most disputed issues in quantum cosmology. Ac- 
cording to Vilenkin's tunnelling proposal ||16||, ^ is given by e +I rather than e -7 . This 



choice of sign is inconsistent with Eq. (1.1), as it leads to the inverse result for the 



pair creation rate: Exp = 1/Tnbp- In our case, we would get Ttp = exp(+7r/A e ff). 
Thus black hole pair creation would be enhanced, rather than suppressed. De Sitter 
space would be catastrophically unstable to the formation of black holes. Since the 
radius of the black holes is given by 1/ y/A C Q, the black holes would be more likely 
the larger they were. Clearly, the tunnelling proposal cannot be maintained. On 
the other hand, Eq. ( |3.15|) , which was obtained from the no boundary proposal, 
is physically very reasonable. It allows topological fluctuations near the Planckian 
regime, but suppresses the formation of large black holes at low energies. Thus the 
consideration of the cosmological pair production of black holes lends strong support 
to the no boundary proposal. 



4 Classical Evolution 

We shall now consider neutral black holes created at any value 0o > 1 of the scalar 
field and analyse the different effects on their evolution. Before we take quantum 
effects into account, we shall display the classical solution for a universe containing 
a pair of black holes. We shall demonstrate explicitly that it behaves according to 
the First Law of black hole mechanics. 
With a rescaled inflaton potential 

Vffl = -^m 2 2 , (4.1) 

87T 

the effective cosmological constant will be 

A eff = m 2 2 . (4.2) 

In the previous section we learned that the black hole radius remains constant, at 
1/VA, in the Lorentzian regime. But this was for the simple model with fixed A. 
The effective cosmological constant in Eq. ( |4.2|) is slightly time dependent. Thus we 
might expect the black hole size to change during inflation. 
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Indeed, for < t < approximate Lorentzian solutions are given by |L8| 



0(f) = <f) -mt, (4.3) 



a(t) = — — cosh 



t 



m / 0(f) df 



o 



(4.4) 



6(0 = — (4-5) 

m(f>{t) 

d s 2 = -dt 2 + a(t) 2 dx 2 + b(t) 2 dtt 2 2 . (4.6) 

Since we are dealing with a degenerate solution, the radii and r c of the black hole 
and cosmo logical horizons are equal: 

r h = r c = b(t). (4.7) 

According Eq. ( f4.5| ) they will expand slowly together during inflation as the scalar 
field rolls down to the minimum of the potential V and the effective cosmological 
constant decreases. At the end of inflation they will be approximately equal to m" 1 . 

One can think of this increase of the horizons as a classical effect, caused by 
a flow of energy-momentum across them. If the scalar field were constant, its 
energy-momentum tensor would act exactly like a cosmological constant. The flow 
of energy-momentum across the horizon would be zero. However, the scalar field 
is not constant but is rolling down hill in the potential to the minimum at <fi = 0. 
This means that there is an energy-momentum flow across the horizon equal to 

M = A T ah l a l b = ^, (4.8) 

where A = Airb 2 is the horizon area, T a ^ is the energy-momentum tensor for the 
massive scalar field, given by 

T ab = ^d a <P d b <P - ^-g ab (d c <j) d c <P + m 2 2 ) , (4.9) 
47r Sir v ' 

and l a is a null vector tangent to the horizon: 

One would expect the horizons to respond to this flow of energy across them by 



an increase in area according to the First Law of black hole mechanics [21 



M — —A, (4.11) 



12 



where k is the surface gravity of the horizon. We will show that this equation is 
indeed satisfied if the horizon growth is given by Eq. (|4.5| ). 

The values of k for general Schwarzschild-de Sitter solutions are derived in the 
Appendix. Due to the slow change of the effective cosmological constant we can 
approximate the surface gravity at any time t in our model by the surface gravity 
in the model with a fixed cosmological constant A = A c fj(t). In the degenerate case 
which we are considering now, k will thus be given by 



Eq. ( [4. 11| ) becomes 



k = JA cS = nuj). (4.12) 



md> d ( 4-7T \ 1 . 

= < 413 > 

which agrees with Eq. ( |4.8| ). 

It should be pointed out that this calculation holds not only for the black hole 
horizon, but also for the cosmological horizon. Moreover, an analogous calculation 
is possible for the cosmological horizon in an ordinary inflationary universe without 
black holes. Thus, in hindsight we understand the slow growth of the cosmological 
horizon during inflation as a manifestation of the First Law of black hole mechanics. 



5 Quantum Evolution 

So far we have been neglecting the quantum properties of the inflationary spacetime 
presented above. It is well known that in a Schwarzschild-de Sitter universe, radi- 



ation is emitted both by the black hole and the cosmological horizon [p2j] . To treat 
this properly, one should include the one-loop effective action of all the low mass 
fields in the metric g^ u . By using a supersymmetric theory one might avoid diver- 
gences in the one-loop term, but it would still be impossibly difficult to calculate in 
any but very simple metrics. Instead, we shall use an approximation in which the 
black hole and cosmological horizons radiate thermally with temperatures 

T b = ^, T c = ^. (5.1) 
2tt' 2tt K J 

This quantum effect must also be included in the calculation of the energy flow 
across the horizons. For the saddle point metric, Eq. ( |4.6| ), it has no consequence: 
in the Nariai solution the black hole and cosmological horizons have the same radius 
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and surface gravity. Thus they radiate at the same rate. That means they will be 
in thermal equilibrium. The black holes will not evaporate, because they will be 
absorbing as much as they radiate. Instead, their evolution will be governed by the 
classical growth described above. 

However, the Nariai metric is an idealization. (Strictly speaking, it does not even 
contain a black hole, but rather two acceleration horizons.) Due to quantum fluctua- 
tions there will be small deviations from the saddle point solution, corresponding to 
a Schwarzschild-de Sitter spacetime which is not quite degenerate. The radius b of 
the two-sphere will not be exactly constant along the S 1 , but will have a maximum 
r c and a minimum rt>, which we identify with the cosmological and the black hole 
horizons, respectively. (The other black hole horizon of the pair will lie beyond the 
cosmological horizon and will not be visible in our universe.) Since the black hole 
horizon is slightly smaller, it will have a higher temperature than the cosmological 
horizon. Therefore the black hole will radiate more than it receives. There will thus 
be a net transfer of energy from the smaller horizon to the larger one. This will 
cause the larger horizon to grow faster, and the smaller one to shrink until the black 
hole vanishes completely. We show below that black holes can still grow with the 
cosmological horizon until the end of inflation, if they are either created sufficiently 
large, or start out very nearly degenerate. However, we shall see that none of these 
conditions is easily satisfied. 

Black holes created during the final stages of inflation will survive until the end 
of inflation simply because they will be relatively large and cold. One can estimate 
the minimum size they must have by treating them as evaporating Schwarzschild 
black holes, which have a lifetime on the order of M 3 , where M is the mass at which 
the black hole is created. In terms of the value of the scalar field at creation, M ~ 
60 = (m0o) -1 - Inflation ends after a time of (fio/m. Therefore black holes created 
at 0o — m~ 1 / 2 will certainly survive until the end of inflation. They would continue 
to grow slowly during the radiation era, until the temperature of the radiation falls 
below that of the black holes. They will then start to evaporate. By Eq. ( p. 15 ), 
however, such black holes will be suppressed by a factor of 



One must therefore investigate the possibility that the small black holes, which 
can be created in abundance, start out so nearly degenerate that they will grow 
with the cosmological horizon until they have reached the "safe" size of m~ 1//2 . We 
need to determine how nearly equal the horizon sizes have to be initially so that the 
black hole survives until the end of inflation. 




(5.2) 
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If we take the thermal radiation into account, the flow across the horizons now 
consists of two parts: the classical term due to the energy flow of the scalar field, as 
well as the net radiation energy transfer, given by Stefan's law. Applying the First 
Law of black hole mechanics to each horizon, we get: 



Kb 
8tt 

H c 

8tt" 



2 2 

m ty 



2 2 

m r„ 



- (aA h T* - aA c T c 
+ (aA b T b 4 - aA c T c 



(5.3) 
(5.4) 



where a = 7r 2 /60 is the Stefan-Boltzmann constant. Using Eq. flS.lp , we obtain two 
coupled differential equations for the horizon radii: 



1 

hZ r T p 



2 2 

m r b — 



240tt 
1 



,2,4 



„2.,4 



r h K h r c n c 



2 2 

m r n + 

c 240vr 



.2 .,4 



,2„4 



r b K b r c K c 



(5.5) 
(5.6) 



The exact functional relation between the surface gravities and the horizon radii 
is generally non-trivial. However, the above evolution equations can be simplified 
if one takes into account that a nucleated black hole pair must be very nearly 
degenerate if it is to survive until the end of inflation. We can therefore write 



r h (t) = b(t)[l-e(t)\, r c (t) = b(t)[l + e (t)\ 



(5.7) 



with e <C 1. The surface gravities can also be expressed in terms of e (see the 
Appendix); to first order they are given by: 



1-3*) 



(5i 



We assume that b{t) behaves as in Eq. ( |4.5| ) for the Nariai solution as long as 
e(t) 1. Eqs. Q5.5D , ( |5.6|) are then identically satisfied to zeroth order in e. In the 
first order, they give an evolution equation for e: 



-m b 



This equation can be integrated to give 



e{t) = e 



2/3 
I ex P 



720tt 



1 1 

180tt63 




(5.9) 



(5.10) 
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For the unsuppressed, Planck size black holes we have 0o — m ~ x ■ If they grow 
with the horizon, they will reach the safe size, which corresponds to = m^ 1 ' 2 , 
after a time 

tsafe = m~ 2 (l - m 1 ' 2 ) . (5.11) 

If a Planck size black hole is to have survived until this time, i.e. if e(t sa f e ) < 1, then 
the initial difference in horizon sizes may not have been larger than 

-2 

(l - m 2 ) 



e - ax = m V3 exp 



m 



720tt 



(5.12) 



The probability P(eo < e™ ax ) that the two horizons start out so nearly equal 
obviously depends on the distribution of the initial sizes of the two horizons. The 
semi-classical treatment of the quantum fluctuations which cause the geometry to 
differ from the degenerate case, for general values of the effective cosmological con- 
stant, is an interesting problem by itself, and beyond the scope of this paper. We 
hope to return to it in a forthcoming paper on complex solutions in quantum cosmol- 
ogy. However, here we are working at the Planck scale, so that the semi-classical 
approximation will break down anyway. It therefore seems reasonable to assume 
that the initial sizes of the horizons are distributed roughly uniformly between zero 
and a few Planck lengths. This means that 

P (e < 6- ax ) « ^ ax « exp i-^A ■ (5-13) 



Since m ~ 10 we conclude by comparison with the suppression factor (|5.2|) that 
it is considerably less efficient to create Planck size black holes which would grow 
to the safe size than just to create the large black holes. Both processes, however, 
are exponentially suppressed. 



6 Charged Black Holes 

Although nearly all neutral black holes will evaporate during inflation, those that 
have a magnetic charge won't be able to because there are no magnetically charged 
particles for them to radiate. They can only evaporate down to the minimum mass 
necessary to support their magnetic charge. Let us therefore introduce a Maxwell 
term in the action and re-examine the pair creation of primordial black holes: 

/ = / fog 1 ' 2 (R-2A- F^ v ) - -L / d\h l ' 2 K. (6.1) 

167T JM, 87T JE 
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The S 2 x S 2 bubbles in spacetime foam can now carry magnetic flux. The action 
of the Maxwell field will reduce their probability with respect to neutral bubbles. 
Thus magnetically charged black holes will also be pair created during inflation, 
though in smaller numbers than neutral black holes. However, once created, they 
can disappear only if they meet a black hole with the opposite charge, which is 
unlikely. So they should still be around today. 

We shall estimate the number of primordial charged black holes present in the 
observable universe. For the purpose of calculating the pair creation rate during in- 
flation, we can use the solutions for a fixed cosmo logical constant, since A e g changes 
slowly. There exists a three-parameter family of Lorentzian charged Schwarzschild- 
de Sitter solutions. They are usually called Reissner-Nordstrom-de Sitter solutions 
and labeled by the charge q and the "mass" \x of the black hole, and by the cosmo- 
logical constant A: 

ds 2 = -U(r)dt 2 + U(r)- l dr 2 + r 2 dQ 2 2 , (6.2) 

where 

U(r) = l - ^ + 4 ~ W- ( 6 - 3 ) 
r r z 6 

We are interested in the cases where the black holes are magnetically, rather than 
electrically charged. Then the Maxwell field is given by 

F = qsinOdO Ad(f). (6.4) 

In an appropriate region of the parameter space, U has three positive roots, which 
we denote, in ascending order, by r- l} r Q and r c and interpret as the inner and outer 
black hole horizons and the cosmological horizon. They can serve as an alternative 
parametrisation of the solutions. For general values of q, /i and A the metric ( p.2|) 
has no regular Euclidean section. The black holes which can be pair created through 



regular instantons lie on three intersecting hypersurfaces in the parameter space []23 
as seen in Fig. 2. They are called the cold, the lukewarm and the charged Nariai 
solutions. In the cold case, the instanton is made regular by setting r; = r D , which 
corresponds to an extremal black hole. The lukewarm hypersurface is characterised 
by the condition q = fi. It corresponds to a non-extremal black hole which has the 
same surface gravity and temperature as the cosmological horizon. This property 
is shared by the charged Nariai solution, which has r Q = r c . For q = 0, its mass is 
given by \x = 1/ (3\/A) and it coincides with the neutral Nariai universe we discussed 
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1/4 



Figure 2: The hypersurfaces in the parameter space of the Reissner-Nordstrom- 
de Sitter solution which admit regular instantons. The plot is of the dimensionless 
quantities fi 2 A vs. q' 2 A. The points N and D represent the uncharged Nariai and 
de Sitter solutions. The curve DC corresponds to the lukewarm instantons. At C 
it meets the curve NU of charged Nariai instantons. The cold instantons lie on the 
curve CU; the point U represents the two ultracold solutions. 



earlier. For larger charge and mass, it is still the direct product of two round two- 
spheres, however with different radii. The cold, lukewarm and de Sitter solutions 
all coincide for q = /i = 0. The largest possible mass and charge for the lukewarm 
solution is q = = 3/(4-\/3A), where it coincides with the charged Nariai solution. 
The largest possible mass and charge for any regular instanton is attained at the 
point where the charged Nariai and the cold hypersurfaces meet. This ultracold 
case has q = 1/(2a/A) and /i = 2/(3v / 2A). It admits two distinct solutions of 
different action. All of these solutions are presented and discussed in detail in the 
comprehensive paper by Mann and Ross |TT[, where the actions are calculated as 
well. We shall now apply these results in the context of inflation. 

Let us consider an inflationary scenario in which the effective cosmological con- 
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stant A e ff = m(f) starts out near the Planck value and then decreases slowly. As in 
the neutral case, one would expect the creation of charged black holes to be least 
suppressed for large A, i.e. at the earliest stage of inflation. However, unlike the 
neutral case a magnetically charged black hole cannot be arbitrarily small since it 
must carry at least one unit of magnetic charge: 

Qo = 7T-, (6.5) 

where eo = \fa. is the unit of electric charge, and a ~ 1/137 is the fine structure 
constant. In the following we shall only consider black holes with q = q , since 
they are they first to be created, and since more highly charged black holes are 
exponentially suppressed relative to them. We see from Fig. 3 that pair creation 




Figure 3: The action of the different types of Reissner-Nordstrom-de Sitter instan- 
tons, as a fraction of the action of the de Sitter instanton, J de sitter = —3ir/A. Note 
that the point U does not correspond to an instanton solution in this diagram; the 
two ultracold instantons are labeled by Ui and U2. E is the instanton with the 
highest pair creation rate during inflation. 

first becomes possible through the ultracold instanton Ui, when A e g has decreased 
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to the value of A u = l/(4g§) 



a:. 



Since this solution has relatively high action, 



however, black hole production becomes more efficient at a slightly later time. It 
will then occur mainly through the charged Nariai instanton, which has lower action 
than the cold black hole. When 



(6.6) 



A cfr = A E = (Ay/3 - 6) a 
the pair creation rate reaches its peak and is given by 



r = exp [- (/ E - /, 



de Sitter 



exp 



2 + v/3 



7T 



2a 



(6.7) 



by Eq. ( |1.1| ). As A eff decreases further, the pair creation rate starts to decrease and 
soon becomes vastly suppressed. For 



< A cfr < A c = 

the lukewarm instanton has the lowest action, 



4 ' 



3tt 



'luke 



A. 



7T4 



off 



aA, 



eff 



(6.8) 



(6.9) 



It will therefore dominate the black hole production until the end of inflation. 

If we ask how many charged black holes were produced during the entire infla- 
tionary era, we have to take into account that the density of the black holes pair 
created at the early stages of inflation will be reduced by the subsequent inflation- 
ary expansion. As a consequence most of the charged primordial black holes in our 
universe were produced near the end of inflation, as we shall show here. The number 
of such black holes per Hubble volume pair created during one Hubble time when 
<p = pc is given by 



exp [- (/luke - J de Sitter)] = exp 



3 7T 



a mcj) 



pc 



(6.10) 



At the end of inflation, these black holes have a number density per Hubble volume 
of 

<4nd(0pc) 



r ( ( / , pc)</'pc eX P 



2#o 



(6.11) 
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where the cubic factor is due to the growth of the Hubble radius and the exponential 
factor reflects the inflationary expansion of space. 

Since we are dealing with exponential suppressions, practically all primordial 
magnetically charged black holes in the universe were created at the value of pc 
which makes d cnd maximal. This occurs for 



imax 
pc 



7T 



1/6 



3 am 2 



(6.12) 



so that the approximate total number per Hubble volume at the end of inflation is 
given by 



D 



end 



iend 



max\ 
pc J 



71 



1/2 



3 am 2 



exp 



9tt 2 



16 am 2 



1/3 



(6.13) 



We take the values of the Hubble radius to be 



H cnd 



m 



H. 



eq 



10 54 , H~ 



10 



59 



(6.14) 



respectively at the end of inflation, at the time of equal radiation and matter density, 
and at the present time. Therefore since the end of inflation the density has been 
reduced by a factor of 10~ 91 m~ 3 / 2 and the Hubble volume has increased by 10 177 m 3 . 
We multiply -D C nd by these factors to obtain the number of primordial black holes 
in the presently observed universe: 



n ~ 



10 



86 



'2 \ V 2 



3 a 



exp 



9tt 2 



1/3- 



16 am 2 



(6.15) 



With a = 1/137 and m = 10~ 6 , the exponent will be on the order of — 10 5 . Thus, 
in ordinary Einstein-Maxwell theory, it is very unlikely that the observed universe 
contains even one magnetic black hole. 

However, in theories with a dilaton both the value of the electric charge and 
the effective Newton's constant can vary with time. If both were much higher in 
the past, the effective value of am 2 would not be so small and the present number 
of magnetic black holes could be much higher. We are currently working on this 
question. 
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7 Summary and Conclusions 



Since the quantum pair creation of black holes can be investigated semi-classically 
using instanton methods, it has been widely used as a theoretical laboratory to 
obtain glimpses at quantum gravity. The inflationary era is the only time when 
we can reasonably expect the effect to have taken place in our own universe. We 
chose to work in a very simple model of chaotic inflation, which allowed us to expose 
quite clearly all the important qualitative features of black hole pair creation. The 
inflationary universe was approximated as a de Sitter solution with a slowly varying 
cosmological constant. Similarly, neutral black holes produced during inflation were 
described by a degenerate Schwarzschild-de Sitter solution. Their pair creation 
rate was estimated from the no boundary proposal, by comparing the probability 
measures assigned to the two solutions. We found that Planck size black holes are 
plentifully produced, but at later stages in inflation, when the black holes would 
be larger, the pair creation rate is exponentially suppressed. This fits in with the 
usual instanton prescriptions for pair creation. The tunnelling proposal, on the other 
hand, fails to make physically reasonable predictions. The consideration of black 
hole pair creation thus lends support to the no boundary proposal. 

We analysed the classical and quantum evolution of neutral primordial black 
holes. Classically, the black hole horizon and the cosmological horizon have the 
same area and temperature. The two horizons grow slowly during inflation. We 
showed that this is due to the classical flow of scalar field energy across them, 
according to the First Law of black hole mechanics. Quantum effects, however, 
prevent the geometry from being perfectly degenerate, causing the black hole to be 
hotter than the cosmological horizon. As a consequence, practically all neutral black 
holes evaporate before the end of inflation. 

Finally, we turned to magnetically charged black holes, which can also be pair 
created during inflation. Even if they have only one unit of charge, they cannot 
evaporate completely and would still exist today. The pair creation rate is highest 
during the early stages of inflation, when the effective cosmological constant is still 
relatively large. Black holes created at that time, however, will be diluted by the 
inflationary expansion. Most of the charged primordial black holes were therefore 
created near the end of the inflationary era, where they would not be diluted as 
strongly. However, in Einstein-Maxwell theory they are so heavily suppressed that 
we must conclude that there are no primordial black holes in the observable universe. 
It will therefore be interesting to examine inflationary models that include a dilaton 
field. One would expect to obtain a much higher number of primordial black holes, 
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which may even allow us to constrain some of these models. 
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Appendix 

In this Appendix we show how to calculate the surface gravities of the two hori- 
zons in the Schwarzschild-de Sitter solution. This space-time possesses a regular 
Euclidean section only in the degenerate (Nariai) case, where the two horizons have 
the same radius. Neutral black holes pair created during inflation will therefore 
start out nearly degenerate. We present a suitable coordinate transformation for 
the nearly degenerate metric, introducing a small parameter e, which parametrises 
the deviation from degeneracy. The surface gravities and Euclidean action are cal- 
culated to second order in e, yielding a negative mode in the action. We explain 
why our results differ from those obtained in Ref. ||15|| . 

The Lorentzian Schwarzschild-de Sitter solution has the metric 

ds 2 = -U{r)dt 2 + U{r)- l dr 2 + r 2 d{l 2 2 , (A.l) 

where 

U(r) = 1 - — - -Ar 2 . (A.2) 
r 3 

For < /i < |A _1//2 , U has two positive roots and r c , corresponding to the 
cosmological and the black hole horizon. The spacetime admits a timelike Killing 
vector field 

K = 4 t , (A.3) 
where jt is a normalisation constant. The surface gravities Kb and k c , given by 



K h c = lim 



{K a V a K b ) (K C V C K 



1/2 



(A.4) 
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depend on the choice of 74. To obtain the correct value for the surface gravity, 
one must normalise the Killing vector in the right way. In the Schwarzschild case 
(A = 0) the natural choice is to have K 2 = — 1 at infinity; this corresponds to 7$ = 1 
for the standard Schwarzschild metric. However, in our case there is no infinity, 
and it would be a mistake to set j t = 1. Instead one needs to find the radius r g 
for which the orbit of the Killing vector coincides with the geodesic going through 
r g at constant angular variables. This is the two-sphere at which the effects of 
the cosmological expansion and the black hole attraction balance out exactly. An 
observer at r g will need no acceleration to stay there, just like an observer at infinity 
in the Schwarzschild case. One must normalise the Killing vector on this "geodesic 
orbit" . Note that this is a general prescription which will also give the correct result 
in the Schwarzschild limit. It is straightforward to show that 



A 



so that 



It = U(r g ) 



-1/2 



Eq. (|A~4D then yields 



2JU(r g ) 



(9A/i 2 ) 
dU 



1/3 



-1/2 



dr 



(A.5) 



(A.6) 



(A.7) 



r=r b 



In order to consider the pair production of black holes, we need to find a Eu- 
clidean instanton which can be analytically continued to the metric ( |A.1| ). The 
obvious ansatz is 

ds 2 = U(r)dr 2 + U{r)- l dr 2 + r 2 d£l 2 2l (A.8) 

where r is Euclidean time. Again one can define a constant 7 T which will normalise 
the timelike Killing vector on the geodesic orbit. In order to avoid a conical singu- 
larity at a horizon one needs to identify r with an appropriate period r ld , which is 
related to the surface gravity on the horizon by 



.id 



27T7 r / K. 



(A.9) 



Usually only one of the two horizons can be made regular in this way, since 
their surface gravities will be different. They will be equal only for \i = |A~ 1//2 , 
when the two roots of U coincide. In this degenerate case one can remove both 
conical singularities simultaneously and obtains a regular instanton. As was first 
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pointed out in [|I5" |, the fact that = r c does not mean that the Euclidean region 
shrinks to zero. The coordinate system (|A.8 ) clearly becomes inappropriate when 
\i approaches its upper limit: the range of r becomes arbitrarily narrow while the 
metric coefficient U{r)~ l grows without bound. One must therefore perform an 
appropriate coordinate transformation. If we write 



9/i 2 A = 1 - 3e 2 , < e < 1, 



(A.10) 



the degenerate case corresponds to e 
coordinates ip and x by 



0. We then define new time and radial 



1 



if>; r 



1 2 4 „ 
1 + ecosx - -e + -e cos% 
o 9 



(A.11) 



With this choice of ip we have 
vector 



K a 



A to second order in e, so that the Killing 




dip 



(A.12) 



has unit length on the geodesic orbit. We have chosen the new radial coordinate 
X so that U vanishes to forth order in e for cosx = ±1. This is necessary since U 
contains no zero and first order terms and we intend to calculate all quantities to 
second non-trivial order in e: 



u(x) 



• i i 
sin x e 



-e cos y H — e cos y H — e 
3 A 3 A 9 



(A.13) 



Thus the black hole horizon corresponds to x = 71 an d the cosmological horizon to 

x = o. 

The new metric obtained from the coordinate transformations (|A.11[) 



is 



ds 2 



1 



2 2 1 

-e cos x + -e 2 cos 2 x - -e 2 ) sin z x dip' 



+ 



1 

A 
1 

A 

1 , . 

— ( 1 + 2e cos x + e cos x — 77 e ' 



1 + -ecosx 



9 



e cos x rfx 



rffi 2 . 



(A. 14) 



In the degenerate case, e = 0, this is the Nariai metric: the topological product of 
two round two-spheres, each of radius 1/yfK. There are no conical singularities if 
the Euclidean time ip is identified with a period 2tt. For general e the two horizons 
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cannot be made regular simultaneously; it is clear from a geometrical standpoint 
that ip must be identified with period 



d_ 

X=0, 7T I 



a/ 9i>i> 



2tt 



x=o> v I 



2 

1 ± -e 

3 



6 



(A.15) 



to prevent a conical singularity at the cosmological or black hole horizon, respec- 
tively. 

One can calculate the surface gravities k c and using the Euclidean version 
of Eq. ( |A.4| ) and the Killing vector ( |A.12j ); equivalently, one could use the relation 
K = 27T7^/^ ld to obtain the same result: 



K c, b 



, 2 11 2 

1 =F -e H e 2 

3 18 



(A.16) 



This equation is useful for the analysis of the radiation energy flux in a nearly 
degenerate Lorentzian Schwarzschild-de Sitter universe, since each horizon radiates 
approximately thermally with the temperature T = k/2tt. 

We will now calculate the Euclidean action of the metric ( |A.14| ) and show that it 
posesses a negative mode in the direction of decreasing black hole mass. The total 
instanton action is given by 



AV A r 5 r A h 5; 



b«b 



57T 



57T 



•>7T 



(A.17) 



where V is the four- volume of the geometry, A Cj b are the horizon areas and 5 Cj b are 
the conical deficit angles at the horizons. Of course, all of these quantities depend on 
the value we choose for ip ld . Obvious options are: to leave it at 2tt even in the non- 
degenerate case, thus introducing a deficit angle on the cosmological horizon and an 
excess angle at the black hole horizon, or as an alternative, to make one of the two 
horizons regular, thereby causing a larger excess or deficit at the other horizon. The 
most interesting of these cases is the one in which we choose a regular cosmological 
horizon. In this case, the metric (|A.14| ) will lie on the interpolation between the 
Euclidean Nariai and de Sitter universes, since the latter has only a cosmological 
horizon, which ought to be regular. The Euclidean actions for these universes are 
— 27r/A for the Nariai, and — 37r/A for the de Sitter. Since no intermediate solution 
is known, one would expect the action to decrease monotonically as one moves away 
from the Nariai solution. In other words, this particular perturbation of the Nariai 
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metric should correspond to a negative mode in the action. Indeed, if we identify ijj 



with the period ip ld , the action in Eq. ( [A. 171) turns out as 



17tt 



A 9A 



+ 0(e 4 



(A.18) 



(The same result is obtained for the period ip^, while for ip ld = 2tt the negative 
mode is given by — 207re 2 /9A.) 

The coordinate transformations, the perturbed Nariai metric and the negative 



mode given here differ from Ref. |y| for various reasons. The authors did not ensure 
that U = on the horizons, and the Killing vector wasn't renormalised properly. 
Also, they identified Euclidean time with period 2n even in the non-degenerate case, 
which is not appropriate to the physical situation we are trying to analyse. 
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